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A search is performed in symmetric 3-jet hadronic Z-decay events for evidence of 
colour singlet production or colour reconnection effects. Asymmetries in the angular 
separation of particles are found to be sensitive indicators of such effects. Upper 
limits on the level of colour singlet production or of colour reconnection effects are 
established for a variety of models. 
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1 Introduction 



The term 'rapidity gaps' denotes regions of angular phase space devoid of particles. They are 
expected in low-pr diffractive processes, where separate colour singlet hadronic systems are 
produced, well separated in phase space and associated with either the target or projectile 
particles. Events containing large rapidity gaps, attributed to colour singlet exchange or colour 
reconnection effects, are also observed, in association with high-p^ jet production, at Hera pQ 
and at the Tevatron As shown in Figures^, and^), by crossing symmetry, similar gaps 
may be expected in three-jet hadronic Z decays at LEP. The corresponding diagram has four 
final state partons, but, because of its generally low energy, the two jets associated with the 
colour singlet object are typically unresolved. 

Large rapidity gaps are observed in ~10% and 1-2% of events with two high-p^ jets at 
Hera and at the Tevatron, respectively. In electron-hadron or hadron-hadron collisions 
particles produced by the spectator partons of the underlying event frequently destroy large 
rapidity gaps associated with the primary hard scattering process. The associated 'gap survival 
probability' is estimated [3] to be about 20% at the Tevatron. An advantage of the Z decay 
study is the absence of this suppression factor, as there is no underlying event. A disadvantage 
is that the maximum possible size of the angular gap is smaller for Z decays compared to ep 
or pp collisions. Particle and energy flow in the inter-jet regions have been studied [I] using 
three-jet events in e + e~ annihilations. These studies revealed that the region between two- 
quark jets have lower particle and energy flows relative to naive expectation from independent 
fragmentation models. This was also observed in studies which compared three-jet events with 
two-jet events having a hard photon in the final state [5]. 

The analysis presented in this Letter is performed with hadronic Z-decay events recorded 
by the L3 detector [Oj using 75.14 pb -1 of data from the 1994-1995 Z-pole data taking periods. 
In the new method presented here, a search is performed for gaps in angular phase space 
in symmetric 3-jet hadronic Z-decay events. The method exploits the different particle flows 
between quark and antiquark jets and either the quark or antiquark jet and the gluon jet. After 
removing particles near to the jet cores, angular gaps between particles in the inter-jet regions 
are analysed, and various asymmetry variables are formed, as detailed below. 

Studies of rapidity gaps in hadronic Z decays using as variable the pseudo-rapidity of par- 
ticles relative to the thrust axis were previously reported A recent study [H] used the axes 
of tagged gluon jets. The analysis presented in this Letter, based on global event variables, is 
complementary to this study in the sense that the jet cores are excluded from the analysis so 
as to minimise fragmentation effects whereas previous analyses rather used rapidity gaps as a 
tool to investigate the details of gluon fragmentation [H]. 

The present analysis extends the notion of comparing particle and energy flow in the region 
opposite to the quark jets as well as in the region between the two-quark jets by introducing 
new asymmetry variables that are sensitive to the relative difference in colour flow between all 
the inter-jet regions. 

A first application of these asymmetries is to exploit differences in colour flow between 
events where colour singlet systems are produced (Csp) and conventional gluon colour octet 
production (Cop). As shown in Figure Hh, in Cop, colour flow is present between the qg and 
qg gaps and is inhibited by destructive interference in the qq gap. As shown in Figure HJl, the 
colour string in the Csp is drawn between the quark and the anti-quark so that an appreciable 
colour flow occurs also in the qq gap. 

A second application is to investigate colour reconnection (CR) effects. Partons originating 
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from a hard scattering process are eventually transformed into hadrons and this hadronisation 
process requires specification of the colour flow pattern among the partons. In the simplest 
models, the colour flow associated with the final state partons is fixed during the hard scattering 
process. However, there may be subsequent rearrangement of the colour flow. At the pertur- 
bative level this requires the exchange of at least two gluons between the partons. Coloured 
strings, normally stretched between a quark and a gluon as shown in Figure can be rear- 
ranged in the colour reconnection picture so as to create colour singlet quark pair in association 
with a colour singlet gluon pair, whose colour strings then hadronise independently, as shown 
in Figure To study CR effects, the Gal [9 model as well as CR as implemented in Ari- 
adne [TUj and Herwig [TTJ are considered. 

Studies of the determination of the W boson mass using fully hadronic W-pair decays, 
indicate CR effects as the dominant source of theoretical systematic uncertainty. If the same 
CR algorithm is valid for both Z and W-pair decays, limits on the level of CR effects, established 
experimentally at the Z-pole, can be used to constrain the systematic uncertainty on the W- 
mass determination. The present analysis is thus complementary to the direct measurement of 
CR effects in hadronic decays of W-pairs [TJ] . 

2 Event and Particle Selection 

Well balanced hadronic Z-decay events are selected by cuts on the number of calorimetric clus- 
ters with energy greater than 100 MeV, N c i uster , on the total energy observed in the calorimeters, 
E vis , the energy imbalance along the beam direction, Eh, and the energy imbalance in the plane 
perpendicular to the beam direction, Ej_. The cut on the number of calorimetric clusters re- 
jects low multiplicity events such as t~t + final states. About two million hadronic Z-decay 
candidates are selected. 

Symmetric three-jet events with a jet-jet angular separation of about 120° are then selected 
using the Jade algorithm [T3], with the jet resolution parameter set to 0.05. The angles between 
jets i and j, 0^, are required to be within ±30° of the symmetric topology. Using the selection 
criteria: 

N duster > 12, 0.6 < E ms /^7s < 1.4, E\\/E vis < 0.40, 

E ± /E vis < 0.40, 0i2,0 2 3,03i e [90°, 150°], 

where \fs is the centre-of-mass energy, about 70000 three-jet events are obtained. In order 
to distinguish quark jets from gluon or colour singlet jets, the energy-ordered quark jets are 
tagged by cuts on the b-tag discriminant 1 >: 

D{ et > 1.25, D 3 2 et > 1.25 D{ et < 1.5 

As shown in Figure these cuts strongly enhance the gluon fraction in jet 3. This selection 
tags 2668 events with a gluon purity of 78%. 

To study the particle flow, calorimetric clusters are selected which satisfy at least one of the 
following criteria: 

• energy greater than 100 MeV in the electromagnetic calorimeter (ECAL) and at least 900 
MeV in the hadron calorimeter (HCAL); 

^The jet b-tag discriminant of jet i, containing n tracks, is defined as: Df et = — log 10 P where P = 
P™ X^j=o ( — m -P™)Vi! an d P" — n?=i-^!- Here, Pj is the probability that the jth track in the jet orig- 
inates at the primary vertex. 
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• energy greater than 100 MeV in the ECAL with a minimum of 2 crystals hit; 

• energy greater than 1800 MeV in the HCAL alone. 

These cuts reject noisy clusters and take into account different thresholds in the calorimeters. 
The distributions of the cluster multiplicity with these selection criteria show good agreement 
between data and Monte Carlo, with residual differences below 2.5%. 

3 Monte Carlo Samples 

The Jetset Parton Shower (PS) Monte Carlo program is used to model Cop. Two simple 
models are used to simulate the expected colour flow in Csp: events of type qq7 are generated 
with a photon effective mass as in the gluon jet mass distribution. In the first model, Csl, the 
photon is replaced by a boosted di-quark jet. In the second model, Cs2, the photon is replaced 
by a gluon fragmenting independently. The total particle multiplicity for both these models 
agree with Jetset within ±1 unit. 

For CR studies, the Gal model, implemented in the Pythia Monte Carlo program |15j . 
uses a default value of 0.1 for the colour recombination parameter, Rq. This value is obtained jU] 
by fitting the model to HI data on the diffractive proton structure function. For this study 
the fragmentation parameters of the model are tuned to Z-decay data 2 ) for three different 
values of the colour recombination parameter: Rq = 0.05, 0.1, 0.2. The Ariadne and Herwig 
generators, with and without CR, are also tuned to Z-decay data to determine their basic 
fragmentation parameters. The colour reconnection probability in Herwig is set to its default 
value of 1/9 jllj . Similarly, the parameters affecting colour reconnection in Ariadne are kept 
at their default values [Ej, para(26)=9 and para(28)=0. 

4 Inter-Jet Gap Asymmetries 

After selection of three-jet events, the particle momenta are projected onto the event plane 
defined by the two most energetic jets. In order to minimise the bias from jet fragmentation, 
particles in a cone of 15° half-angle about the jet axis direction are excluded from the analysis. 
The angles of the remaining particles are measured in this plane with respect to the most 
energetic jet. In order to achieve uniformity in the event-to-event comparison, these angles 
are rescaled so as to align jets at 0°, 120° and 240°. This is achieved by scaling the angle of 
a particle to its nearest jet by the ratio between 120° and the opening angle of the two jets 
between which the particle is located. 

Two gap angle definitions are used [TKj . as shown in Figures andEb: the minimum angle, 
Bij, of a particle measured from the bisector in the gap ij, and the maximum separation angle, 
Sij, between adjacent particles in the gap. In Figures Ek and Eli the minimum energy of the 

^The QCD models are tuned using several global event shape distributions at ~ mz' the minor on the 
narrow side the jet resolution parameter for the transition from 2- to 3-jet in the Jade ^2] algorithm, the 
fourth Fox- Wolfram moment ^3 an d the charged particle multiplicity. For models implemented in Pythia, 
the tuned parameters are the QCD cut-off parameter, A, and the string fragmentation parameters, b and 
<tq, affecting longitudinal and transverse components of the hadron momenta. In Herwig, the QCD cut-off 
parameter and the parameters controlling hadronisation, CLMAX (maximum cluster mass) and CLPOW (the 
power of the mass in the expression for the cluster splitting criterion) are tuned. 
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calorimetric clusters used to define the bisector angle or the maximum separation angle is com- 
pared with the Jetset prediction. Good agreement is obtained, showing that the contribution 
of soft particles, to which rapidity gap distributions are particularly sensitive, is well simulated. 
The angular asymmetry in gap 12 from the Bij angles is defined as 

12 B12 + -B23 + £31 
Asymmetries are also defined from the SV,- angles as 

j^S —S\2 + 5*23 + 531 

S\2 + 5*23 + S31 

The other gap asymmetries: Af-, Af-, with ij = 23, 31 are defined in a similar way. 

Reduced colour flow and thereby larger separation for Csp in gaps 23, 31 with respect to 
gap 12, should thus make Af 2 and Af 2 peak more strongly at positive values for Csp than for 
Cop. The 23 and 31 gap asymmetries of each event are averaged to yield the 'qg' asymmetries 
shown below. 

The angular asymmetry distributions are corrected for detector effects and initial and final 
state photon radiation using bin-by-bin correction factors obtained from events generated with 
the Jetset Parton Shower Monte Carlo program and processed through L3 detector simula- 
tion [IHI- The bin sizes are chosen sufficiently large that migration effects are negligible. The 
correction factors are defined as the ratio of generated particle-level distributions, considering 
all charged and neutral particles, without energy cuts, to the same distributions after detec- 
tor simulation. The particle-level distributions take into account the gluon jet identification 
probability and have a quark flavour composition corresponding to y/s ~ mz- These bin- by-bin 
correction factors for the angular asymmetries typically lie in the range of ±20%. 

The particle-level angular asymmetry distributions of the selected symmetric three-jet events, 
normalised to unit area, are compared to different models, in Figures EJ El and El 

Fractional bin-by-bin systematic uncertainties are estimated by repeating the analysis us- 
ing clusters obtained by combining calorimetric clusters with tracks, as used, for instance, in 
Reference 1121 A variation between 2% and 5% is observed. Furthermore, the cuts on the b-tag 
discriminant are changed so that the gluon purity varies by ±10%, which results in system- 
atic uncertainties between 3% and 8%. Finally, the residual 2.5% difference between data and 
Monte Carlo discussed above is included. These uncertainties are added in quadrature and are 
summarised in Table d The systematic uncertainty due to a change in the jet cone angle cut 
from ±15° to ±20° is found to be negligible. 





Systematic uncertainties (%) 


Variable 


Detector 


B-tag 


Monte Carlo 


Total 


A 12 


5.2 


4.8 


2.5 


7.5 


A qg 


5.9 


3.2 


2.5 


7.1 


A 12 


6.6 


8.1 


2.5 


10.8 


AS 
^qg 


2.8 


4.0 


2.5 


5.5 



Table 1: Systematic uncertainties on the measurements of the asymmetry variables 
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5 Limits on Colour Singlet Production 



As shown in Figure IH the data are in good agreement with the Cop model. The high discrim- 
ination power of the angular asymmetries between the Cop and Csp models is also evident. 
A quantitative comparison is given in Table El The Cop model is in good agreement with the 
data. The Csp models are clearly excluded. As a cross-check, the analysis is repeated using the 
Durham k± algorithm [20J with y cut = 0.01 and y cu t = 0.02 instead of the Jade algorithm, 
effectively defining an independent set of asymmetries. No significant changes are observed. 





x 2 


for 






Variable 


Cop 


Csl 


Cs2 


d.o.f. 




6.4 (0.99) 


356 


262 


19 


A qg 


15.9 (0.60) 


238 


189 


18 




4.8 (0.94) 


1081 


524 


11 


AS 
^qg 


6.7 (0.88) 


334 


266 


12 



Table 2: Values of x 2 obtained from the comparison of the data distributions to colour octet and 
colour singlet models. For Cop the corresponding confidence levels are given in parentheses. 
For the Csp models all confidence levels are less than 10~ 30 . The \ 2 values include systematic 
uncertainties. 



The asymmetry distributions are fitted to a combination of Cop and Csp contributions. 
This is done by minimising a x 2 function defined as 



X\r) = £ 



[fdata r /cSP ^ r ^COP ] 2 



{ a ltat) 2 + { a \yst) 2 



where /' is the content of the zth bin of the distribution, r is the fraction of the Csp component 
and contributions from both statistical, cr\ tat , and systematic, &\ yst) uncertainties are included. 

Good fits are obtained for all the asymmetry distributions. For the variable Af 2 , which 
Monte Carlo studies show to be the most sensitive one, the fit gives r = 0.015 ± 0.024 (stat.) 
± 0.018 (syst.) with x 2 /d.o.f. = 4.5/11 for the Csl model r = 0.025 ± 0.031 (stat.) ± 0.029 
(syst.) with x 2 /d.o.f. = 4.4/11 for the Cs2 model. All the fits give a fraction of events due to 
Csp consistent with zero. The fits to the distributions are then used to obtain a 95% confidence 
level (CL) upper bound on the fraction of Csp events. The asymmetry variables Af 2 and A^ g 
are independent, as are Af 2 and A^ g . These pairs of variables are thus combined in the fits. 
Upper bounds of 6.7% and 10.2% for the Csl and Cs2 models, respectively, are found using 
Af 2 and A^ g . Using Af 2 and A^ g yields slightly weaker limits. 



6 Limits on Colour Reconnection Effects 

The particle-level angular asymmetry distributions are compared to the predictions of several 
different Monte Carlo models in Figures 03 and HI These include the 'no CR' models Jetset, 
Ariadne and Herwig as well as the Gal model and the CR versions of Ariadne and 
Herwig. The x 2 confidence levels (CLs) given by the comparison of these models to the data 
are presented in Table H3 Both the default Gal and Ariadne CR models are excluded by the 
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Af 2 distribution, with CLs of ~ 10 -8 and ~ respectively. The Gal model is also excluded 
by the A^ g distribution which gives CL ~ 10~ 4 . The same distribution has a low CL of ~ 10 -2 
for the Ariadne CR model. However Ariadne, without CR, gives a satisfactory description 
of all of the distributions. Both versions of Herwig are completely excluded, with a best CL of 
~ 10 -8 for no CR and of ~ 10~ 9 for CR, among all of the asymmetry distributions considered, 
suggesting that it cannot be used, with confidence, to simulate the soft hadronisation effects 
that are important for CR studies 3 ^. Consistent results are obtained by repeating the analysis 
using the Durham k± algorithm [2H] with y cut = 0.01 and y cut = 0.02 instead of the Jade 
algorithm. 





No CR 


CR 


Variable 


Jetset 


Ariadne 


Herwig 


Gal 


Ariadne 


Herwig 


A B 


0.99 


0.93 


10- 9 


0.04 


0.27 


io- 11 


A B 
qg 


0.60 


0.13 


io- 8 


10" 4 


0.02 


io- 8 




0.94 


0.80 


1Q -24 


IO" 8 


10~ 6 


IO" 30 


A s 


0.88 


0.78 


io- 8 


0.03 


0.07 


io- 11 



Table 3: x 2 confidence levels obtained from the comparison of the data distributions to differ- 
ent models with and without CR. Low confidence levels are rounded to the nearest order of 
magnitude. 

Fits are performed to the asymmetry distributions to obtain the best value of Rq by inter- 
polating the Monte Carlo distributions with different values of Ro- Good fits are obtained, in 
all cases, with values of Ro consistent with zero. Further fits are then performed to obtain an 
upper limit on _R . Combining the pair of variables Af 2 and A^ g or Af 2 and A^ g , a 95% CL 
upper limit for R of 0.024 is obtained. 

7 Summary and Conclusions 

New observables based on angular separations of particles in the inter-jet regions of symmetric 
three-jet events are introduced and are found to be very sensitive to Csp and to CR effects. 

Upper limits at 95% CL on Csp according to the Csl and Cs2 models of 6.7% and 10.2%, 
respectively, are obtained. Since the fraction of Csp expected on the basis of the Tevatron 
measurements is only 5-10%, after allowing for the effect of gap survival probability, the present 
analysis is not sufficiently sensitive to confirm or exclude a similar effect in hadronic Z decays. 

The Gal model, with the default CR probability, and the Ariadne CR model are unable 
to describe the data. Both the no CR and CR versions of Herwig are completely excluded 
by the data. However, a good description is provided both by Jetset and the no CR version 
of Ariadne. This suggests that the angular asymmetries are also very sensitive to the non- 
perturbative hadronisation model used. Both Jetset and Ariadne have similar, string-like, 
hadronisation models, whereas Herwig uses cluster fragmentation. 

The results presented in this Letter provide important information concerning the systematic 
uncertainty on the W mass resulting from CR effects as estimated by the Gal, Ariadne and 
Herwig Monte Carlo models. For the default value of the Gal CR parameter, Ro = 0.1, the 

3 )This conclusion does not depend on the size of the systematic uncertainties which, even doubled, would still 
give Herwig CLs of ~ 10~ 7 for the Af 2 and A^ g asymmetries. 
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W mass measured from decays of W pairs into four jets is estimated [0] to be shifted by about 
65 MeV. The 0.024 95% CL upper limit obtained in this analysis implies a mass shift of only 
a few MeV. Since the default CR models in Ariadne and Herwig are unable to correctly 
describe the Z-decay data, it is difficult to have confidence in their use to describe CR effects 
in W-pair production. 
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Figure 1: Schematic diagram of colour exchange. Colour singlet propagators are indicated by 
double lines in a) pp and in b) e + e~ reactions. The e + e~ diagram is derived by crossing the 
incoming quark line in the pp diagram. Colour flow is shown by dashed lines for c) Cop and d) 
Csp in 3-jet events from e + e~ annihilation and also e) without and f) with colour reconnection. 
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Figure 3: Definition of a) minimum angle relative to the gap bisector and b) the maximum 
separation angle between adjacent particles for the case of four particles in the sensitive region of 
gap 12. Only particles outside the ±15° cones around the jet axes are considered. Distributions 
of the minimum energy of clusters used to define c) the bisector angle and d) the maximum 
separation in selected symmetric three-jet events compared to the Jetset PS prediction. 
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Figure 4: a) and b) minimum bisector angle gap asymmetries, and c) and d) maximum sepa- 
ration gap asymmetries for gaps 12 and qg, respectively, compared to colour singlet and colour 
octet models. 
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Figure 5: a) and b) minimum bisector angle gap asymmetries, and c) and d) maximum sep- 
aration gap asymmetries for gaps 12 and qg, respectively, compared to models without CR 
effects. 
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Figure 6: a) and b) minimum bisector angle gap asymmetries, and c) and d) maximum sepa- 
ration gap asymmetries for gaps 12 and qg, respectively, compared to models with CR effects. 
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